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ABSTRACT 

We present a ~ 1" (100 pc) resolution 12 CO (3-2) map of the nearby intermediate stage interacting 
galaxy pair NGC 4038/9 (the Antennae galaxies) obtained with the Submillimeter Array. We find 
that half the CO (3-2) emission originates in the overlap region where most of the tidally induced 
star formation had been previously found in shorter wavelength images, with the rest being centered 
on each of the nuclei. The gross distribution is consistent with lower resolution single dish images, 
but we show for the first time the detailed distribution of the warm and dense molecular gas across 
this galaxy pair at resolutions comparable to the size of a typical giant molecular complex. While 
we find that 58% (33/57) of the spatially resolved Giant Molecular Associations (GMAs; a few x 
100 pc) are located in the overlap region, only < 30% spatially coincides with the optically detected 
star clusters, suggesting that the bulk of the CO (3-2) emission traces the regions with very recent 
or near future star formation activity. The spatial distribution of the CO (3-2)/CO (1-0) integrated 
brightness temperature ratios mainly range between 0.3 and 0.8, which suggests that on average the 
CO (3-2) line in the Antennae is not completely thermalized and similar to the average values of 
nearby spirals. A higher ratio is seen in both nuclei and the southern complexes in the overlap region. 
Higher radiation field associated with intense star formation can account for the nucleus of NGC 4038 
and the overlap region, but the nuclear region of NGC 4039 show relatively little star formation or 
AGN activities and cannot be easily explained. We show kinematical evidence that the high line ratio 
in NGC 4039 is possibly caused by gas inflow into the counter-rotating central disk. 
Subject headings: galaxies: individual (NGC 4038, NGC 4039) — galaxies: interactions — galaxies: 
ISM — radio lines: galaxies 



1. INTRODUCTION 

Gravitational interaction between galaxies is now con- 
sidered to be one of the primary mechanisms of galaxy 
evolution. Major mergers, composed of two or more 
galaxies with comparable masses, are widely believed to 
pro vide a way to form elliptical and spheroidal galax- 
ies (jToomre lll977t ). Furthermore, galaxy interaction is 
known to increase star formation activities as seen in the 
increasing fraction of tidally distorted morphologies in 
ultraluminous infrared galaxies (ULIRGs: Lfir > 10 12 
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L Q , SFR - 10 2 M vr- 1 : ISanders fc Mirabell Il996h . In 
addition, the number of merg ing galaxies appear s to in- 
crease at high redshifts (e.g. iBundv et al.l 1200 9!) . The 
obvious first step in characterizing this process and the 
response of the cold gas and its relation to merger in- 
duced starburst activity is to observe nearby merging 
galaxies. 

The Antennae galaxies (NGC 4038/9) is the nearest 
colliding galaxy pair observable from th e northern hemi- 
spher e located at a distance of 22 Mpc (ISchwcizc r et al.1 
2008). Numerical simu lations of the Antennae (e.g. 
iToomre fc Toomira Il972f ) suggested that t he Antennae 
are at the intermediate stage of the collision. Mihos et al.l 
(1993) suggest that the two spiral galaxies passed the ini- 
tial pericenter 210 Myr ago, and the final coalescence will 
take place in about 100 Myr. 

The Antennae have been obser ved at all availabl e wave - 
lcngths from radio to X-ray. iWhitmore et al.l (1999) 
observed the Antennae with the Wide Field Planetary 
Camera 2 on the Hubble Space Telescope (HST) and iden- 
tified thousands of stellar clusters which have formed 
during the merging process. They found that most of 
the star formation occurs in the form of massive compact 
star clusters, which are formed from supergiant molecu- 
lar clouds. Previous observations have found that most of 
the active star forming regions are concentrated in the re- 
gion where the two galaxies overlap. The Herschel-PACS 
maps at 70, 100 and 160 (im revealed that both nuclei 
arc brighter than the HII regions in the arms, but the 
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brightest e mission in three b ands comes from the over- 
lap region (|Klaas et all 12010ft . The mid-infrared (MIR) 
and far-infrared (FIR) emission traces buried star forma- 
tion, which is obscured by dust and therefore is invisible 
at optical wavelengths. They estimated that the integral 
SFR is ~ 22 M Q yr" 1 and the SFR in the southernmost 
knot of the overlap region with an area of 4.68 kpc 2 is ~ 
3 M yr -1 . They also found that several emission knots 
in the overlap region are in different evolutionary stages. 
Low angular resolution observations have revealed strong 
[CII] emission, which arises mainly from photo dissocia- 
tion regions created by far-ultraviol et photons from ho t 
young stars, near the overlap region ([Nikola et al.l [l998). 
In addition, the southernmost knot of the overlap region 
is currentl y the most active area w ith a very young stellar 
complex. Ne ff fc Ulvestadl (|2000ft observed 4 and 6 cm 
radio continuum emission using the VLA and estimated 
that the overall star formation rate seems to be a factor 
of ~ 20 higher than normal galaxies. They found that 
the compact radio source with the strongest thermal ra- 
dio emission requires ~ 5000 05 stars to account for the 
free-free emission. They estimate a supernova rate of ~ 
0.2 yr -1 if the sources with steep radio spectra are dom- 
inated by non-thermal emission produced in supernova 
remnants. 

The distribution and the characteristics of molecular 
gas have been investigated i n detail us i ng mm /submm 
telescopes around the world. iGao et all (|2001ft obtained 
the CO (1-0) map using the NRAO 12 m single-dish tele- 
scope and found that the distribution of the CO (1-0) 
emission is similar to those of the MIR a nd FI R emis- 
sion at a kiloparsec scale. IWilson et al.l (|2000l ) identi- 
fied five supergiant molecular complexes (SGMCs) with 
masses of (3-6) x 10 s M in the overlap region using the 
12 CO (1-0) map obtai ned at the Owen s Valley Radio 
Observatory (OVRO). iZhu et al.l (|2003| ) have obtained 
single-dish maps at the Nobeyama 45m telescope in the 
CO (1-0) line and using the James Clark Maxwel Tele- 
scope (JCMT) in the CO (2-1) and CO (3-2) lines. Each 
nucleus contains a molecular mass of ~ 10 9 M Q and the 
overlap region accounts for approximately 42% of total 
12 CO (1-0) flux of the inner disk of the Antennae. They 
found that the CO integrated intensity ratios are high 
in the overlap region. This is the first case where high 
CO (3-2)/CO (1-0) ratio s are found outsid e a galac- 
tic nucleus. Furthermore, iSchulz et al.l (2007) obtained 
CO (1-0) and CO (2-1) maps using the IRAM 30 m 
Millimeter Radio Telescope and CO (3-2) map using the 
Heinrich Hertz Telescope. The total molecular gas mass 
of the system (~ 10 10 Mq) is about twice the total gas 
mass of the Milky Way. However the molecular gas mass 
within the 1 kpc nucleus of NGC 4038 and NGC 4039 
exceeds that of the central 1 kpc of the Milky Way by 
a factor of almost 100, revealing high gas concentration 
into the two nucleus. ISchulz et al.l (|2007f l also derived 
the line ratios ranging between 0.5 and 0.66 , which are 
lower than the ratios estimated bv lZhu et al.l (|2003l ). 

The purpose of our work is to investigate the physical 
properties and the kinematics of the CO (3-2) emitting 
molecular complexes in the Antennae using ~ 1" (1" 
corresponds to about 107 pc) resolution CO (3-2) obser- 
vations. The CO (1-0) rotational transition has a critical 
density of n cr it ~ 10 2 ' 5 cm 3 and an upper energy level 



E u ~ 5 K, whereas the CO (3-2) rotational transition 
has n crit ~ 10 4 cm 3 and E u ~ 33 K. Thus the CO (3-2) 
emission is a better tracer of denser and warmer molec- 
ular gas than CO (1-0) emission. 

This paper is organized as follows. We describe our 
observations in §2 and our results in §3. In §4, we provide 
a discussion of our results. We present the properties 
and the distribution of identified molecular complexes 
(§4.1), the CO (3-2)/(l-0) brightness temperature ratio 
(§4.2), the kinematics in NGC 4039 (§4.3), and a possible 
molecular gas bubble in NGC 4038 (§4.4). We summarize 
and conclude this paper in §5. 

2. OBSERVATIONS 

We conducted CO (3-2) observations of main bodies of 
the NGC 4038/ 9 with the Submillimeter ArrajQ (SMA; 
IHo et all 12551 on March 2005 and March 2008 in the 
compact configuration and December 2008 in the ex- 
tended configuration. The data were obtained using a 
five-pointing mosaic where the phase centers were cho- 
sen to cover the main CO (3-2) em ission peaks seen in 
the JCMT map of IZhu et all (|2003l ). The coordinates of 
the phase centers are shown in Table Q] and the fields of 
view are shown in Figure Q] (right). The total on source 
observing time was approximately 10 hours. The SIS 
receivers were tuned to the frequency of the CO (3-2) 
line (^ rC st=345.796 GHz) rcdshifted to the systemic ve- 
locity of 1634 km/s. The primary beam of the array at 
this frequency is 32". The correlator had a 2 GHz total 
bandwidth with a 0.8 MHz (0.7 km s _1 ) frequency res- 
olution. The quasar 3C 279 was observed for phase and 
amplitude calibration. Absolute flux calibration was per- 
formed using Uranus. The uncertainty of flux calibration 
is 20 %. 

Data inspection and calibration was carried out using 
the MIR package written in IDL and imaging was done 
using the MIRIAD package. The synthesized beam size 
is l'/42 x 1'.'12, which corresponds to 150 pc x 120 pc 
with a position angle of 29.4°. Natural weighting of the 
visibilities is used to maximize the sensitivity, and the 
achieved rms noise level in 10 km s _1 channel maps is 38 
mjy beam -1 . The continuum emission is not detected 
with an upper limit of 17.7 mjy (3 a). 

3. RESULTS 

3.1. CO (3-2) Distribution and Total Flux 

The CO (3-2) integrated intensity map overlaid on op- 
tical/MLR images is shown in Figure [TJ In addition, we 
compare the CO (3-2) map with multiwavelength images 
from radio and to X-ray (see Appendix A) . In general, the 
locations of the CO (3-2) emitting molecular complexes 
ar e consistent w i th the positions of the SGMCs detected 
bv lWilson etaLl ff2000h in CO (1-0), but the higher angu- 
lar resolution achieved in the CO (3-2) map shows details 
that were unseen in the CO (1-0) image. Both galaxies 
show central CO (3-2) concentrations near the nuclei and 
have one gas arm roughly distributed along the dust lanes 
toward southwest seen in the HST 435 nm image (Figure 
Q](left)). The strongest CO (3-2) emission is associated 

12 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 
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with the nucleus of NGC 4038 with a peak flux density 
(10 km s -1 velocity resolution) of 2.0 ± 0.4 Jy beam -1 . 
The CO (3-2) emission near the nucleus of NGC 4039 is 
composed of at least two distinct large scale components 
separated by 560 pc. The overlap region consists of nu- 
merous emission peaks, with complex structure seen in 
larger (~ 1 kpc) molecular complexes. Spatial correlation 
with optical emission is poor in most of the molecular 
complexes in the overlap region, mostly because these re- 
gions are strongly affected by dust extinction. However, 
the general association with the Spitzer 8 /xm emission 
(Figure Q] (right)), which traces the dusty regions associ- 
ated with CO (3-2) emitting molecular clouds, is much 
better than the optical association. 

The total integrated CO (3-2) flux of main body of 
the Antennae is 4.4 x 10 3 Jy km s -1 . The two nuclei of 
NGC 4038 and NGC 4039 contain 1.4 x 10 3 and 9.4 x 
10 2 Jy km s , respectively. The total flux in the overlap 
region is 2.1 x 10 3 Jy km s , which is approximately 
half of the total integrated flux over the entire galaxy. 
The percentage of missing flux derived by comparing 
the SMA and JCMT flux (jZhu et all [2001 are 55% in 
NGC 4038, 35% in NGC 4039 and 67% in the overlap 
region. The low missing flux in NGC 4039 suggests the 
compact nature of the CO (3-2) emitting clouds, while 
the relatively large missing flux suggests the presence of a 
significant extended (> 3 kpc) component in the overlap 
region. 

The central ~ 100 pc in NGC 4038 is unresolved by 
the ~ 1" beam, but the strong signal allows us to use 
the extended configuration data alone for a higher angu- 
lar resolution image (Figure [2]) . The synthesized beam 
size is 0'.'82 x 0'.'56 by adopting natural weighting of the 
visibilities. Most of the diffuse and extended emission 
below 15(7 in the ~ 1" resolution map made using all 
data (Figure [TJ is resolved out in the high resolution 
map. The deconvolved size of the main component is 
230 pc (2.1") x 150 pc (1.4"), and the FWHM velocity 
width is 110 km s . The brightness temperature Tg 
(the peak intensity) scales with the beam filling factor 
/ and the molecular gas temperature T. From the peak 
intensity of 13 K and assuming that the typical molec- 
ular g as temperature T o f the CO (3-2) emission is < 
30 K (jSchulz et all 120071 ). the beam filling factor is < 
0.4. Thus the distribution of the CO (3-2) emission is 
likely more compact than the region observed with the 
~ 75 pc beam. 

3.2. CO (3-2) Kinematics 

The velocity fields of NGC 4038, NGC 4039, and the 
overlap region are shown in Figure [3] (a), (b) and (c), re- 
spectively. At this spatial resol ution, the velocity a t the 
K s -band nucleus of NGC 4038 (jMengel et al.ll2002h is ~ 
80 km s -1 lower than the surrounding medium, which 
may suggest a rotation or gas streaming into the nuclear 
region along the arm. In contrast, the CO (3-2) emitting 
gas near the nucleus of NGC 4039 shows a large veloc- 
ity gradient (310 km s -1 kpc -1 ) from the north-west to 
south-east, with the velocity ranging from 1560 km s _1 
to 1640 km s _1 . Furthermore, the gas component dis- 
tributed along the optical dust lane of NGC 4039 also 
shows a steep velocity gradient of 170 km s -1 kpc -1 . 
The implications of these kinematic features will be dis- 



cussed in §4.3. In the overlap region, the velocity of the 
gas ranges from Vlsr = 1300 to 1650 km s -1 . The large 
scale kinematics in the overlap region show a north-south 
velocity gradient, which is consistent with t he velocity 
field seen in the lower resolution JCMT map (|Zhu et al.1 
2003). The new SMA map reveals that each molecular 
component shows complex kinematics, with velocity gra- 
dients in excess of 100 km s _1 kpc -1 in each molecular 
complex, some as large as 300 km s -1 kpc -1 . 

The channel maps of the high resolution cube of 
NGC 4038 with 10 km/s velocity resolution are shown 
in Figure|4j The rms noise level (la) is 70 mJy beam -1 . 
At this angular and velocity resolution, the morphol- 
ogy of the emission feature is highly complex around 
the nucleus, with the emission distributed in a semi- 
circle around the known nonthermal radio source domi- 
nated by synchrotron r adiation from supernova remnants 
(|Neff fc Ulvestadl[2000l) . The radio source does not have 
detectable CO (3-2) associated at an rms of 70 mJy 
beam - 1 which corresponds to a molecular gas mass limit 
of 9 x 10 5 Mq. Further discussion on the possible origin 
of this CO (3-2) arc is provided in §4.4. 

4. DISCUSSION 
4.1. Properties of the CO (3-2) Molecular Complexes 

We identified molecular complexes by applying the 
automatic clump id entification algorithm Clumpfind 
([Williams et al.lll994h to the CO (3-2) data. One of pro- 
grams in Clumpfind, CLFIND, works on 3D (RA, Dec, 
velocity) data cubes, searches for local peaks of emis- 
sion, and follows them down to lower intensity levels. 
Then it decomposes the data cube into clumps in which 
the emission is concentrated. While the 2a threshold is 
the reco mmended value for id entifying robust molecular 
clumps (|Williams et al.lfl994T) , our experiments find that 
the total number of identified molecular complexes varies 
by 15% for a 10% deviation from the 2a threshold. This 
is likely because the program falsely detects (~ 20%) 
sidelobes of the real emission. We find that the number 
of identified clumps only changes by a few percent for a 
threshold of 2.2 < a < 2.6. Thus we chose 2.6a for ro- 
bustly detecting molecular clumps with high significance 
in the following analysis. 

A total of 57 molecular complexes are identified, 33 
(58%) of which are located in the overlap region. The 
properties of the molecular complexes such as position, 
systemic velocity (V), radiusv(r), and velocity disper- 
sion (q„) were measu red using the program CLST ATS 
((Williams et al.|["l994fl and are summarized in Table [BTlI 
(Appendix B). We estimate the error in the position, 
radius, and velocity to be < O'.'l, 44 pc, and 5 km s -1 , 
respectively. The radii of the molecular complexes ranges 
between 85 and 348 pc, whereas the velocity dispersion 
ranges between 6 and 36 km s . 

4.1.1. Molecular Mass 
The molecular gas mass is derived by; 

^=3.25xl0 7 X co 

Mq 
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where Xco is the CO luminosity-to-H2 mass conversion 
factor, ScoAv is the velocity integrated flux, v i, s is the 
observing frequency, and Dj, is the luminosity distance 
([Solomon fc Vanden Bout! 1200^. We used X co = 4.8 
M Q /K km s _1 pc 2 ()Solomon fc Barrett! H99ll ) scaled by 
the average CO (3-2) /(1-0) line intensity ratio for each 
region (see §4.2). The molecular masses are estimated 
to be 10 7 - 10 s M©. In addition, the virial mass of the 
molecular complexes are estimated from M V i r = 2rv 2 /G, 
where r and v is the radius and the velocity dispersion of 
the molecular complex, and G is the gravitational con- 
stant. The molecular complex is gravitationally bound 
when their M V i r /Mco is smaller than unity. We find 
that 84 % of the molecular complexes are gravitation- 
ally bound, including the most massive molecular com- 
plex located north of the NGC 4038 nucleus. We note 
that the velocity dispersions used here are likely overes- 
timated by a factor of 2 - 3 (§4.1.3), and hence the de- 
rived virial masses are also significantly overestimated. 
This suggests that almost all of the molecular complexes 
observed at this spatial resolution have M V i r /Mco < 1 
and that they are gravitationally bound. 

4.1.2. Radius and Velocity Dispersion 

The distribution of the radii of the molecular com- 
plexes in each region are shown in Figure [5] The 
Kolmogorov-Smirnov test gave P-value = 0.002 (Tabled]) 
for the comparison of the distribution between NGC 4038 
and NGC 4039, suggesting that the two populations are 
significantly different. There are molecular complexes 
with large radii (> 180 pc ) in NGC 4038, but not in NGC 
4039. The median radii in NGC 4038 and NGC 4039 
are 210 and 120 pc, respectively, and thus the molecular 
complexes near the nucleus of NGC 4038 are on average 
1.8 ± 0.7 times larger than those in NGC 4039. The size 
distribution in the overlap region range between 90 and 
350 pc. The average size of the molecular complexes in 
the southern part of the overlap region (declination < 
-18° 53' 00") is 1.3 times larger than the average size of 
those in the overlap region. In addition, the P-value = 
0.975 and 0.017 (Table [2]) for the comparison of the dis- 
tribution between the southern part of the overlap region 
and NGC 4038, and NGC 4039, respectively. This sug- 
gests that the radii between NGC 4038 an d the south- 
ern pa rt of the overlap region are similar. iKlaas et al.l 
(2010) estimated the SFRs of seven individual emission 
knots from their Herschel maps. Knot Kl, which cor- 
responds to the southern part of the overlap region, has 
the highest SFR (2.78 M Q yr" 1 ) and knot NN, which 
corresponds to the central region of NGC 4038, has the 
second highest SFR (1.17 M Q yr _1 ). Therefore the larger 
molecular complexes, which are possibly formed by clus- 
tering around local peaks of the gravitational potential 
wells, may be used as the signposts of current dusty star 
forming regions. 

The velocity dispersions of the molecular complexes 
in each region are shown in Figure O We conducted a 
Kolmogorov-Smirnov test on the distribution of the ve- 
locity dispersions, and obtained P-value = 0.664 - 0.973 
(Tabled, suggesting that the velocity dispersions among 
the three different regions are similar. In NGC 4039, 
molecular complexes with larger velocity dispersions (25 
- 35 km s _1 ) are located within 500 pc from the nucleus 



and those with a smaller velocity dispersion (10 - 25 km 
s _1 ) are located in the gas arm. Comparing the molec- 
ular complexes in the Antennae with clumps in M64 
which is a molecule-rich sp iral galaxy located at 4.1 Mpc 
(iRosolowskv fc Blitzl [20051) and the Galaxy (|Oka et al.l 
l2007f) . the molecular complexes with the same velocity 
dispersion have an order of magnitude larger molecular 
mass in the Antennae. 

4.1.3. The Relation between Radius and Velocity Dispersion 

On the GMC scales, it is known that the radius i s 
proportional to the vel ocity dispersion ([Larson! [1981|). 
Rosolowsky et al. (2003) find a correlation between the 
radii and velocity dispersions of GMCs in M 33 which is 
a normal spiral galaxy located at 850 kpc, a nd the corre- 
lation appears to hold even i n the outer di sk (Bigic l et al.1 
120101 ). In other instances, iFukui et all (|2008l ) showed 
that the GMCs in the Large Magellanic Cloud (LMC) fol- 
low a similar correlation, but it is offset from that foun d 
for the GMCs in the inner Galaxy ([Solomon et al.lll987l) . 
In order to investigate if such correlation exists in the 
Antennae at GMA scales, we plot the relation between 
the radii and the velocity dispersions of the identified 
molecular complexes (Figure [7]) • The correlation coeffi- 
cients derived from the least squares fit to the different 
regions are 0.68 in NGC 4038, 0.25 in NGC 4039, and 
0.01 in the overlap region, suggesting the absence of a 
strong correlation in each region of the Antennae. At 
these scales, the velocity dispersion likely reflects the rel- 
ative bulk motion of a collection of CO (3-2) emitting 
clumps, and not the intrinsic velocity dispersions per- 
taining to the individual GMCs. Hence the absence of 
a strong correlation at size scales of GMAs is somewhat 
expected, especially in the dynamically disturbed parts 
of the galaxy such as the overlap region. 

We investigate further the relation between radius and 
velocity dispersion of molecular complexes in NGC 4038 
because the correlation coefficient in NGC 4038 is the 
highest among three regions. By using Clumpfind on 
the high resolution data, we identified 11 molecular com- 
plexes that are mostly located 100 pc north of NGC 4038. 
Their radii range between 40 and 80 pc and are about 
the same size as the GMCs in M 33. The relatively tight 
correlation between radius and velocity dispersion is seen 
in the GMCs de tected in a quiescent spiral galaxy M 33 
(|Onoderal[2009[ ) (Figure [7]). However, the velocity dis- 
persions of the molecular complexes in NGC 4038 are 2 
-3 times larger than the GMCs in M 33 and distributed 
between 10 and 30 km s _1 . This is further evidence sug- 
gesting that the velocity dispersion does not reflect the 
intrinsic velocity dispersion but the relative bulk motion 
of CO (3-2) emitting clumps. 

Finally, we note that the lack of significant difference in 
the velocity dispersion of molecular complexes between 
NGC 4038 and NGC 4039 may be a result of limited 
sensitivity and angular resolution. The clump identifica- 
tion routine is unable to separate the apparently smooth 
and continuous distribution of molecular gas into smaller 
clumps. 

4.1.4. Distribution of the Molecular Complexes and Their 
Relation to Star Clusters 

The green circles in Figure [8] show the positions of the 
identified molecular complexes and the crosses represent 
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the positions of 50 most luminous (at V-band) and/or 
the 50 mos t NIR-bright clusters whose ages are older 
than 1 Myr (|Whitmore et alj|2010f >. Overall, only 17 out 
of 57 molecular complexes (30 %) host one or more of 
these star clusters with a median age of 4 Myr estimate d 
from the ages of the star clusters (jWhitmore et al.ll2010h . 
not considerin g projection e f fects and extinction. In a 
related study, IMengel et al.l (|2005fl identified star clus- 
ters using the VLT K s -band and found that about 70 % 
of the star clusters with masses < 10 5 Mq are younger 
than 10 Myr and most of the young clusters with the 
age < 6 Myr a re loc ated in the overlap region. Although 
IMengel et ail (|2005l ) suggest that the CO (1-0) obtained 
at a spatial resolution of ~ 4" in the overlap region shows 
spatial correlation with the location of the star clusters, 
our CO (3-2) emission map with a spatial resolution of 
~ 1" shows a clear offset from the star clusters, assuming 
that the star clusters identified using the V-band image 
overlap with those identified using the K s -band image. 
The offset is evident in the southern part of the over- 
lap region, where the CO (3-2) emission is strong and 
the CO line ratio is relatively high. Hence, the absence 
of a significant spatial correlation between the locations 
of the molecular complexes and star clusters may sug- 
gest that the molecular gas traced in CO (3-2) emission 
corresponds to the regions of future star formation, or re- 
gions where star formation occurred in the past ~ 1 Myr 
because no star cluster which is younger than 1 Myr cor- 
respond to the molecular complexes. This timescale is 
similar to the typical timescale for disass ociation from 
the m olecular cloud which is under 3 Myr (|Lada fc Ladal 
l200l . 

4.2. Distribution of the integrated brightness 
temperature ratio 

We investigated the integrated brightness temperature 
CO (3-2)/(l-0) ratios using the SMA CO (3-2) and 
OVRO CO (1-0) data ([Wilson et alJl200Ch . While we 
assume the same beam filling factor for both transitions, 
it is possible that the true distribution of the CO (3-2) 
emitting cloud is more clumpy and concentrated than 
the clouds emitting CO (1-0) emission. Therefore, the 
peak line ratios we estimate here are lower limits. We 
clipped the visibilities so that both data have the same 
shortest uv range (minimum uv distance = 10.48 kA = 
9.13 m) and then convolved the CO (3-2) data to the an- 
gular resolution (5".49 x 3".84) of the CO (1-0) image 
before estimating the line ratios. Finally, the line ratios 
are clipped at the 3<r level. 

The resultant ratio map is shown in Figure HJ We 
separate the ratio map into five regions (NGC 4038, 
NGC 4039, Complex 1, Complex 2, and Complex 3) 
and estimate the mean ratio for each region. The mean 
integrated intensity ratio ranges from 0.3 to 0.6 (Ta- 
ble [3]). Here the error bars are derived from the un- 
certainties in the flux calibration. On average, the CO 
line in the Antennae is not thermalized up to J = 3-2. 
Th ey are slig h tly le ss than the average ratios estimated 
by iZhu et al.l j2003) using single dish telescopes, which 
ranged between 0.67 and 0.91. In contrast, the line ratios 
are close to un ity in the souther n part of Complex 3 and 
the two nuclei. Ilono et all (|2009l ) derived 0.93 for the cen- 
tral regions of 12 dusty U/LIRGs in the local universe, 
which is also in agreement with the starburst galaxies 



Arp220 and M82 where the derived ratios are > 1 (e.g. 
iGreve et al.ll2009f) . Thus, the average excitation condi- 
tions analyzed at ~ 1 kpc scale in the Antennae are close 
to the average properties of nearby galaxies ( e.g. 0.44 is 
the a verage ratio among 7 nearby galaxies; iMao et al~l 
l2010f ). except for two localized regions where the excita- 
tion conditions are similar to the central regions of star- 
burst galaxies. 

Finally, we mention two issues that likely affect the 
analysis and interpretation of the line ratios. One is- 
sue is the inapplicability of the Rayleigh- Jeans ap- 
proximation, especially for the higher frequency CO (3- 
2) line. In such case the brightness temperature ra- 
tio can be le ss than unity (0.8 - 0.9 for a temperature 
of 40 K; see IHarris et al.l §010)), assuming a thermal- 
ized and optically thick, single component cloud. The 
other possibility arises from the filling factor difference 
between the CO (3-2) and CO (1-0) lines. Such a multi- 
component cloud may be more realistic than adopting 
a single component cloud, and higher angular resolution 
ALMA/ AC A observations of both lines should provide 
us with the exact distribution of each line tracers in the 
future. 

4.2.1. The Nucleus of NGC 4038/9 

The mean 12 CO (3-2)/(l-0) ratio of NGC 4038 is 0.6 
± 0.2, which is the highest among the averages of the 
five complexes. The spatial distribution of the ratio pre- 
sented in Figure [5] shows that the ratio is 0.75 in the 
center, and this is two times higher than the ratios in 
the outer region. The regions with particularly high line 
ratios correspond well with the peaks of the MIR and 
NIR emission, qualitatively suggesting the correlation 
between high line ratio and embedded star formation. 
This is consistent with observations in nearby galaxies 
where the ratio appears t o increase from the outer disk to 
the central region s (e.g. Ilsraelll2009l IWilson et all 120091 : 
iBoone et al.l 1201 if ). A similar trend is also seen in the 
Milky Way galaxy, where the ratio increases from ~ 0.5 
in the inner disk to ~ 0.9 towards th e galactic cente r 
where the star formation is more active (|Oka et al.f 2007). 
Furthermore, signatures of more recent star formation 
activity in NG C 4038 are seen in opt i cal and X-ray ob- 
servations (e.g. iWeedman et al.l 2005; Zc zas et alJ 12006b 
IKlaas et alJl20ia T suggesting the coexistence of a vari- 
ety of star clusters formed at various stages of the tidal 
activity. 

The line ratio near the nuclear region of NGC 4039 
is high (~ 1), despite a significantly lower mean ratio 
(0.5±0.1). The star formation activity in the central re- 
gion of NGC 4039 is only modest (SFR = 0.50 M yr" 1 ; 
IKlaas et al.ll2010l ). which is only half of the value derived 
in NGC 4038. The absence of the Br 7 line in the stel- 
lar continuum is an indication of the dominance of old 
giants and red supe rgiants, further suggesting inactivity 
(Gilbert et all 120001 ) . Thus it is not likely that the high 
line ratio in the nuclear region is produced by current 
star formation. 

The CO line could be collisionally heated and excited 
up to high energy levels if NGC 4039 ha rbors an AGN 
(jMatsushita et al.112004 iKrips et alJl201lD . A bright dif- 
fuse X-ray emission region is associated with the K s -band 
nucleus within 220 pc, and also with the region with high 
line ratio. The radio continuum has a steep spectrum, 
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which suggests non-thermal emission a rising from super- 
nova remnants (INeff fc Ulve stadl [2000l). or X-rav binaries 
(jZezas et al.ll2002allb[ ). Furthermore, TBrandl et ail (|20Q9h 
found that the observed infrared continuum is consis- 
tent with a pure starburst model and concluded that the 
spectrum of the nucleus is produced by emission from 
dust that is heated by star formation alone. Thus, these 
multi-wavelength properties of NGC 4039 provide evi- 
dence against the presence of an AGN. 

Another possible cause for the high line ratio is gas 
collision arising from violent gas kinematics. As shown 
in §3, the distribution of the CO (3-2) gas in NGC 4039 
is composed of two main components; the nuclear con- 
centration and the single-arm gas that appears to stream 
along the optical dust lane. Although the spatial reso- 
lution of the ratio image is limited to 500 pc, the re- 
gion where the two components connect is consistent 
with the position that displays the highest ratio. The 
properties of the molecular complexes located at this re- 
gion are markedly different from the average properties 
of NGC 4039, where the radii of the two molecular com- 
plexes are larger. In addition, the velocity dispersion of 
the molecular complexes in this regions are ~ 50% larger 
than the average of NGC 4039 . Furthermore, the reces- 
sion velocities (see Table IB.1[) of these molecular com- 
plexes are consistent with the velocity of the clouds with 
the highest ratio. This observational evidence suggests 
a scenario where the molecular gas flows into the central 
disk (see §4.4), collides, and results in high excitation. 
Future high resolution observations of shock tracers is 
necessary to confirm this scenario. 

4.2.2. The Overlap Region 

The mean 12 CO (3-2)/(l-0) ratio is 0.4± 0.1 in Com- 
plex 1, 0.3± 0.1 in Complex 2, and 0.5± 0.1 in Complex 
3. Although the average ratios are all consistent within 
the errors, gradients in the line ratio are seen in each 
complex (Fi gure El), with a parti cularly steep gradient in 
Complex 3. IZhang et all (|2010l ) derived the FUV - 24 
/im broad-band spectral energy distribution (SED) and 
found that the UV/optical emission gradually become 
redder and weaker from Complex 1 (which corresponds 
to their region 6) to Complex 3 (region 4). This suggests 
that the ratio of young star population (< 10 Myr) to 
intermediate/old star population (> 10 Myr) increases 
from Complex 1 to 3. The continuum emi ssion at MIR 
wavele ngths in Complex 1 (SGMC 1 in iWilson et al.l 
(2000)) and Complex 2 (SGMC 2) are flat, but Com- 
plex 3 (S GMC 3-5) contains s trong and rising continuum 
emission (Wil son et al.l [2000h . The mid- infrared contin- 
uum emission in Complex 3 can be produced by massive 
O stars. Similarly, strong FIR continuum emission seen 
in Complex 3 (Kl) has the h ighest SFR and S FE and 
maximum dust temperature (jKlaas et al.1 l2010i ). This 
observational evidence suggests that Complex 3 is likely 
to be a very young complex and star forming histories 
are different from Complex 1 and 2. This is consistent 
with the higher CO line ratio in Complex 3. 

The overlap region of the Antennae was observed with 
a sing le pointing in the 12 CO (2-1) line at the SMA. 
From our analysis of the archival data, the distribution 
of the CO (3-2) and CO (2-1) maps are similar overall, 
but subtle differences are seen particularly in Complex 3. 
We found that the CO (2-1) peak in Complex 3 is located 



100 pc north of the CO (3-2) peak, and the CO (1-0) 
peak is further displa ced toward th e nort h by 200 pc 
as shown in Figure [TU] Bro gan et al.l (|2010ft detected an 
H2O maser within 100 pc of the CO (3-2) peak, and esti- 
mated from the VLA 3.6 cm continuum that the ionized 
gas is equivalent to 2000 - 5000 O* stars. Furthermore, 
the most ma ssive cluster is embedd ed 350 pc southeast of 
Complex 3 (|Whitmore et al.ll2010D . These are observa- 
tional indicators that suggest the earliest stages of active 
star formation in the southern part of Complex 3 are 
entangled with more mature clusters that were formed 
1 Myr ago. The high line ratio coincides with this re- 
gion, suggesting that heating of the interstellar medium 
due to UV emission from newly born stars is responsible 
for collisionally exciting the CO gas to the J=3-2 level. 

4.3. Evidence for a Counter Rotating Nuclear Disk and 
Gas Inflow Toward NGC 4039 

The spatial distribution of the CO (3-2) emission in 
NGC 4039 (Figure [11] (left)) clearly shows two distinct 
molecular gas components that are separated by a pro- 
jected distance of 560 pc. One component (component 
A; see Figure UTTleft)) is apparently distributed around 
the nucleus and the other component (component B; see 
[TIT left)) is distributed along the dust lane seen in the B- 
band image. The masses of components A and B are 5.2 
x 10 s and 4.5 x 10 8 M Q , respectively. The PV diagram 
of component A is shown in Figure [TT] (right). The sys- 
temic recession velocity of NGC 4039 determined from 
our data is 1610 km s . The mean gas velocity of com- 
ponent A is blue- and redshifted by -50 and +30 km 
s -1 with respect to the K s -band nucleus to the west and 
east side, respectively. Then assuming that component A 
is in the plane of the galaxy, this suggests that compo- 
nent A rotates around the galactic center of NGC 4039. 
We compare the dynamical mass and molecular mass of 
component A and find that the disk is dynamically sta- 
ble, assuming that component A is a disk with a rota- 
tion velo city of 40 km s" 1 and t h e disk inclination of 
i = 60 ° (jToomre fc Toomrdll972l ). E owne s & Solo mon! 
(1998) found a high occurrence of rotating 100 pc scale 
nuclear disks or rings from their molecular gas obser- 
vations of ULIRGs. While the kinematic signature of 
a disk in NGC 4039 is less obvious than the ULIRGs 
in lDownes fe; Solomon! (|1998l ) mainly because of a lower 
S/N ratio, presence of such nuclear disks may be a com- 
mon phenomenon in the nuclear regions of colliding sys- 
tems. A clear characteristic seen in NGC 4039 is that 
this component appears to rotate in the opposite direc- 
tion from the galactic rotation defined by the direction 
of the spiral arms seen in optical images, which is defined 
from the direction of the trailing arm. 

iV-body/smoothed p article hydrodynamic (SPH) sim- 
ulation performed by Barnes & Hcrnquist (1996) pre- 
dicts that a non-axisymetric potential can enhance the 
stellar torques exerted on the gas, causing gas inflows into 
the central regions in a merging process. These inflows 
may be kinem atically dist i nct, a s shown in the simula- 
tion studies bv llono et al.l (|2004l ). Component B shows 
the predicted kinematic signature of gas streaming along 
the arm (~ 1 kpc), as shown in Figure [3] In addition, 
a high line ratio is seen at the location where Compo- 
nents A and B connect and the velocity dispersion of 
the molecular complexes are ~ 50% larger than the av- 
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erage of NGC 4039 as shown in §4.2.1. If Component B 
is indeed related to gas inflow, it is then possible that 
high gas excitation is caused by collisions between the 
inflowing gas component and the counter-rotating disk. 

4.4. A possible molecular gas bubble in NGC 4038 

The high resolution image (~ 0"7; Figure [2]) allows us 
a detailed investigation of the molecular gas surround- 
ing the nuclear region of NGC 4038. The molecular and 
virial masses of the central component are 3.2 x 10 8 M Q 
and 2.0 x 10 8 M , respectively, suggesting that this com- 
ponent is gravitationally bound. Furthermore, by us- 
ing Clumpfind on the high resolution data, we identified 
11 molecular complexes that are mostly located 100 pc 
north of the galaxy center. The virial ratios for all of 
the clumps are less than unity, suggesting gravitation- 
ally bound clouds. 

A nonthermal radio source with an energy equivalent 
to 115 times Cassiopeia A is located 110 pc north of 
the K s -band peak (square sign in Figure [5]) (Neff & Ul- 
vestad 2000). The distribution of CO (3-2) appears to 
form an arc centered around this radio source, with an 
apparent absence of significant CO (3-2) emission at 
the radio source. Similar molecular bubb les have been 
found in M 82 (e g. Matsushita et al]|200l and NGC 253 
(|Sakamoto et al.H2006[) . Thev attributed it to supernova 
blasts that swept away the molecular gas from the sur- 
rounding medium. Following iSakamoto et al.l ()2006f ). we 
estimated the kinetic energy needed to create the bubble 
using, E » (107r/3)pi; 2 -R 3 , where p = lAnomn {no is 
the hydrogen number density of the prebubble medium, 
and rnu is the the atomic mass of hydrogen) is the initial 
mass density of the interstellar matter, R is the size of a 
bubble, and v is the expansion velocity of molecular gas. 
We derived R from the distance between the nonthermal 
radio source and the CO (3-2) peak in Figure[5J which is 
110 pc, and used the offset from systemic velocity at the 
CO (3-2) peak as v of 30 km s _1 . The estimated energy 
is 9 x 10 46 J, adopting no = 10 2 cm~ 3 , which is a typical 
density of interstellar matter in the Galaxy and nearby 
galaxies. This energy corresponds to the kinetic energy 
released by 10 - 100 supernova explosions and also to 
the released energy estimated from the radio continuum. 
Thus a high frequency of recent supernova explosions is 
a possible cause of the apparent absence in the CO (3-2) 
distribution. 

5. SUMMARY 

We present interferometric 12 CO (J = 3-2) obser- 
vations toward main bodies of the Antennae galaxies 
(NGC 4038/9) obtained using the SMA. The image 
shows gas concentration at both nuclei and the overlap 
region, and the emission is roughly distributed along the 
dust lanes seen in the B-band image. The integrated in- 
tensity over main body of the Antennae is 4.4 x 10 3 Jy 
km s -1 and the overlap region contains almost half of the 
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emission. We detect using the Clumpfind algorithm 57 
molecular complexes, 33 of which located in the overlap 
region. Adopting a galactic CO-to-H2 conversion factor 
and the estimated CO (3-2)/(l-0) line intensity ratios, 
the molecular complexes appear to have masses of 10 7 - 
10 s M Q . While we find that 58 % (33/57) of the spatially 
resolved GMAs are located in the overlap region, only 
< 30% are spatially coincident with the star clusters de- 
tected in the optical and/or NIR images, suggesting that 
the bulk of the CO (3-2) emission traces the regions with 
very recent or near future star formation activity. 

We find that the GMAs in the Antennae do not follow 
the well known radius- velocity dispersion relation obeyed 
by GMCs in Milky Way and other nearby galaxies. The 
molecular complexes with large radii are distributed in 
NGC 4038 and the southern part of the overlap region, 
where the star formation rates are high. Therefore, the 
larger molecular complexes, which are possibly formed 
by clustering around local peaks of the gravitational po- 
tential wells, may be used as the signposts of current 
dusty star forming regions. In NGC 4039, the molecular 
complexes with the larger velocity dispersion (25-35 km 
s _1 ) are located within 500 pc of the nucleus and those 
with the smaller velocity dispersion (10-25 km s _1 ) are 
located in the gas arm. 

We estimated the CO (3-2)/(l-0) brightness tempera- 
ture ratio using the SMA CO (3-2) and OVRO CO (1-0) 
data, assuming a same beam filling factor for both transi- 
tions. The mean ratios range between 0.3 and 0.6, which 
suggests that, on average, the CO line in the Antennae is 
not thermalized up to J = 3-2. However the line ratios 
in two nuclei and the southern of the overlap region are 
close to unity in some localized regions. The high ratio 
in the nucleus of NGC 4038 and the overlap region is 
likely caused by star formation because the line ratio be- 
comes high for clouds that are close to active star forming 
regions. In contrast, the nucleus of NGC 4039 has a rel- 
atively inactive star forming region. One possibility for 
causing the high line ratio in the nucleus of NGC 4039 
is a presence of a hidden AGN. The other is gas colli- 
sion between the inflowing gas and the counter-rotating 
central disk. 

The Atacama Large Millimeter/Submillimeter Array 
(ALMA) will begin scientific observations in the second 
half of 2011. The Antennae is one of the best targets 
observed using the ALMA because it passes near zenith 
at the ALMA site and the beam pattern will be a nearly 
circular shape. Using the superb sensitivity and angu- 
lar resolutions offered by ALMA, we can observe other 
molecular lines including high- J CO in the Antennae and 
study the properties of molecular gas quantitatively. 

We thank the anonymous referee for useful comments. 
We also thank Sachiko Onodera for providing the data 
of GMCs in M 33 and giving helpful comments, which 
improved our discussion. 
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Fig. 1. — The CO (3-2) integrated intensity (between 1304 and 1804 km s x ) contour map overlaid on the HST 435 nm image (left) 
and the Spitzer 8 (im image (right). The contour levels are 6 Jy km s _1 X 2, 3, 5, 7, 9, 11, 13, 15, 20, 25. The galaxy and region names are 
labeled in the left figure. The circles in the right figure show the primary beam centered on each pointing. The ellipse in the lower right 
corner of both figures shows the synthesized beam size (l"42 X 1'.'12). These integrated intensity maps are clipped at 2<r level using the 
MPS task, MOMNT. 



TABLE 1 

The phase centers 



RA (J2000) Dec (J2000) 



FOV 1 


12:01:53.45 


-18 


53:12.40 


FOV 2 


12:01:54.61 


-18 


53:08.00 


FOV 3 


12:01:55.46 


-18 


52:55.98 


FOV 4 


12:01:55.05 


-18 


52:43.96 


FOV 5 


12:01:53.01 


-18 


52:01.68 



TABLE 2 

KS-TEST P-VALUE 



regions 


radius velocity dispersion 


(NGC 4038, NGC 4039) 


0.002 


0.771 


(NGC 4038, overlap region) 


0.199 


0.664 


(NGC 4039, overlap region) 


0.006 


0.973 


(NGC 4038, the southern part of the overlap region) 


0.975 


0.675 


(NGC 4039, the southern part of the overlap region) 


0.017 


0.997 
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Fig. 2. — CO (3-2) integrated intensity contour map in NGC 4038. The black contours show the CO (3-2) emission using only the 
data observed in the extended configuration. The contour levels are 5 Jy km s — 1 X 5, 7, 9, 11. The small ellipse in the lower right 
corner indicates the higher angular resolution (0''82 X 0''65). The grey contour map is the CO (3—2) map made with all the data. The 
contour levels are 6 Jy km s _1 X 5, 7, 9, 11, 13, 15, 20, 25. The large ellipse in the lower right corner indicates the angular resolution 
(1.42" X 1.12 "). The cross sign sho ws the galactic center defined by the VLT K s -band image and the square sign shows the nonthermal 
radio source INeff fc UlvestadlboOOT) . 



TABLE 3 

The mean integrated intensity ratio of each region 



mean ratio 


NGC 4038 


0.6 ± 0.2 


NGC 4039 


0.5 ± 0.1 


Complex 1 


0.4± 0.1 


Complex 2 


0.3± 0.1 


Complex 3 


0.5± 0.1 
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FlG. 3. — The velocity field of NGC 4038 (upper-left), NGC 4039 (upper-right), and the overlap region (bottom). The cross signs in the 
two upper figures show the galactic center defined by the VLT K s -band image. 
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Fig. 4. — The channel maps in NGC 4038 using only the data observed in the extended configuration. The contour levels are noise level 
(Ict = 70 mjy beam ) X 3, 4, 5, 6, 7, 8, 9, 10. The cross sign shows the galactic center and the box sign shows the radio source. 
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Fig. 5. — Distribution of the radii of the detected molecular complexes in NGC 4038 (left), NGC 4039 (middle) and overlap region 
(right). 
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Fig. 6. — Distribution of the velocity dispersions of the detected molecular complexes in NGC 4038 (left), NGC 4039 (middle) and 
overlap region (right). 
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Fig. 7. — The relation betwee n the radius and the velocity dispersion of the identified molecular complexes. The data of clumps in M33 
was obtained bvlOnodcra (200j3). 
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Fig. 8. — Distribution of the molecular complexes and their relation to star clusters. The green circle signs show the position of the 
identified molecular complexes. The circle sizes represent the size of the molecular complexes. The blue, red and orange cross signs show 
the location of star cl usters classified into the 50 most luminous (at V-band), most IR-bright clusters and both category, respectively 
(Whitmorc ct al. 2010). The five white circle show the primary beam centered on each pointing. 
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Fig , 9. — The integrated brightness temperature CO (3-2)/(l-0) ratios. The background image is the HST 435 nm (Whitmore et al. 

rj993). 
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Fig. 10. — The closeup of the Complex 3. The Spitzer 8 lira contour map overlaid on the HST Ha image. The contour levels are 40, 
50, 60, 70, 80, and 90 % of the peak intensity. The upper, middle, and lower circles sh ow the CO (1-0), th e CO (2-1), and the CO (3-2) 
peaks, respectively. The upper cross shows the lo cation of the H2O mas er observed by Brogan ct al. (2010) and the lower cross shows the 
location of the most massive cluster identified by Whitmorc ct al. (2010). 
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Fig. 11. — The right figure is the PV diagram along the green cut line shown in the left figure. The position whose corresponds to the 
VLT Ks-band peak of NGC 4039 shown the green circle in the left figure. The contour levels are 6 Jy km a -1 X 2, 3, 5, 7, 9. 
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APPENDIX 
A. MULTIWAVELENGTH IMAGES 

We compare the CO (3-2) map with multiwavelength images from radio and to X-ray. Figure lA~T1 shows the CO (3- 
2) contour maps of five regions (NGC 4038, NGC 4039, Complex 1, Complex 2, and Complex 3) overlaid on 14 different 
images. 

The CO (1-0) map was taken with the OVRO interferometer bv lWilson et al.l (|2000f ). The stron gest emissio n in al l 
regions is located in NGC 4038 and the emission in Complex 2 is strongest in the overlap region, felaas et al.l ([2010) 
conducted FIR observations at 70, 100 and 160 /im using the Herschel-PACS. The brightest emission in three bands 
comes from Complex 3. Complex 1 and Complex 2 together become the brightest area in the 160 /x m map. The 
24 fjm and four near -infrared (NIR) a t 3.6, 4.5, 5.8, and 8.0 /im images were obtained with the MIPS (|Rieke et al.l 
2004) and the IRAC ()Fazio et al.ll2004l ) on board the Spitzer Space Telescope, respectively. We use the basic calibrated 
data obtained using the LEOPARD software. The strongest emission in all bands, especially 24 /im, is located in 
Complex 3. The Ha i mage ( F656N), and the B-band image (F435W) have been observed using the ACS on board the 
HST (jWhitmore et all feOlO). We use the archival data obtained from the MAST (Multimission Archive at STScI). 
There is no optically bright emission in Complex 2 and Complex 3. Both the far-ultraviolet (FUV; ~ 1516 A) and 
the near-ultraviolet (NUV; ~ 2267 A ) images were taken with the Galaxy Evolution Explorer (GALEX) Ultraviolet 
Atlas of Nearby Galaxies distributed bv lGil de Paz et aTl J2007). The strong UV emission is detected west of the peak 
in Complex 3, where the br i ght o ptical, NIR and MIR emission is also observed. We obtained the archival X-ray 
image which iFabbiano et al.1 (j2003| ) observed with Chandra. The bright X-ray emission is associated with the nucleus 
of NGC 4039. 
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Fig. A.l. — The back ground images are the CO (1-0), the 160^tm, the lOOfim, the 70^tm, the 24/im, the 8.0/^m, and the 5.6^tm images 
from left to right. The contour maps are CO (3—2) emission and the contour levels are 6 Jy X 2, 3, 5, 10, 15, 20, 25. The red cross signs 
show the CO (3—2) peaks each region. The length of each side corresponds to 1.6 kpc. Note that the coordinates of the 160/^m, the 100/im, 
and the 70 fim maps have a few arcseconds deviations from the coordinates of the CO (3-2) map. 
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Fig. A.l. — Continued. The background images are the 4.5/xm, the 3.8^tm, the Ha, the B-band, the NUV, the FUV, and the X-ray 
images from left to right. 
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B. CATALOG OF THE MOLECULAR COMPLEXES 

TABLE B.l — The molecular complexes 

TD RA □ 2000) Dec (.12000) V r " 5^ — WZ, W. 



